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A METHOD FOR PREDICTING FULL SCALE BUFFET RESPONSE 


WITH RIGID WIND TUNNEL MODEL FLUCTUATING PRESSURE DATA 

VOLUME I 

PREDICTION METHOD DEVELOPMENT AND ASSESSMENT 


Sy 

Atlee M. Cunningham, Jr. , David B. Benepe, 
Darlene Watts, and Paul G. Waner 
Fort Worth Division of General Dynamics 


Abstract 


This report dociaments the development and assessment of a 
method with which fluctuating pressure data obtained from rigid 
scaled wind-tunnel models can be used to predict flexible full- 
scale buffet response. The method requires unsteady aerodynamic 
forces, natural airplane modes, and the measured pressure data 
as input. A gust response computer program is used to calculate 
buffet response due to the forcing function posed by the measured 
pressure data. By calculating both S3rmmetric and antisymmetric 
solutions, upper and lower bounds on full-scale buffet response 
are formed. Final results are given in the form of upper and 
lower bounds on the power spectral densities and the RMS values 
for angle of attack variation in maneuvers at several Mach- 
altitudes . Comparisons of predictions with flight test results 
are made and the effects of horizontal tail loads and static 
aeroelasticity are shown. Discussions are also presented on 
the effects of primary wing torsion modes, chordwise and span- 
wise phase angles, and altitude. 

This first volume presents a description of the prediction 
method, how the method was evolved, and comparisons of predicted 
rms responses with flight test data for assessing the capability 
of the method. Included with the comparisons are complete 
descriptions of the flight test maneuver and data points, air- 
plane conditions, and calculated modes and wind-tunnel data 
points used in the prediction. The second volume presents the 
power spectral density plots made for all of the predictions. 


Abstract (continued) 


The supporting analysis of F-lllA flight data was conducted 
in two phases. Phase I concentrated on wing and fuselage buffet- 
ing response at a wing leading edge sweep of 26 degrees and sub- 
sonic speeds and supported the initial development of the the 
prediction method. That work is reported in three contractor 
reports . 

Benepe, D. B. , Cunningham, A. M. , Jr. , and Dunmyer , W. D. : 

An investigation of Wing Buffeting Response at Subsonic 
and Transonic Speeds: Phase I F-lllA Flight Data Analysis. 

Volume I - Summary of Technical Approach, Results 
and Conclusions, NASA CR- 152109. 

Volume II - Plotted Power Spectra, NASA CR- 152110. 

Voltame III - Tabulated Power Spectra, NASA CR- 152111. 

Phase II extended the analysis effort to wing leading edge 
sweeps of 50 and 72.5 degrees at both subsonic and low supersonic 
speeds and provided data for assessing horizontal tail response 
and evaluating the final prediction method. The Phase II effort 
is reported in three additional contractor reports. 

Benepe, D. B. , Cunningham, A. M. , Jr., Traylor, S. , Jr., 
and Dunmyer, W. D. : An Investigation of Wing Buffeting 

Response at Subsonic and Transonic Speeds : Phase II 

F-llIA Flight Data Analysis. 

Volume I - Summary of Technical Approach, Results 
and Conclusions, NASA CR- 152112. 

Volume II - Plotted Power Spectra, NASA CR- 152113. 

Volume III - Tabulated Power Spectra, NASA CR- 152114. 
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A METHOD FOR PREDICTING FULL SCALE BUFFET RESPONSE 


WITH RIGID WIND TUNNEL MODEL FLUCTUATING PRESSURE DATA 

VOLUME I 

PREDICTION METHOD DEVELOPMENT AND ASSESSMENT 


By 

Atlee M. Cunningham, Jr. , David B. Benepe, 
Darlene Watts, and Paul G. Waner 
Fort Worth Division of General Dynamics 


Abstract 


This report documents the development and assessment of a 
method with which fluctuating pressure data obtained from rigid 
scaled wind-tunnel models can be used to predict flexible full- 
scale buffet response. The method requires unsteady aerodynamic 
forces, natural airplane modes, and the measured pressure data 
as input. A gust response computer program is used to calculate 
buffet response due to the forcing function posed by the measured 
pressure data. By calculating both symmetric and antisymmetric 
solutions, upper and lower bounds on full-scale buffet response 
are formed. Final results are given in the form of upper and 
lower bounds on the power spectral densities and the RMS values 
for angle of attack variation in maneuvers at several Mach- 
altitudes. Comparisons of predictions with flight test results 
are made and the effects of horizontal tail loads and static 
aeroelasticity are shown. Discussions are also presented on 
the effects of primary wing torsion modes, chordwise and span- 
wise phase angles, and altitude. 

This first volume presents a description of the prediction 
method, how the method was evolved, and comparisons of predicted 
rms responses with flight test data for assessing the capability 
of the method. Included with the comparisons are complete 
descriptions of the flight test maneuver and data points , air- 
plane conditions, and calculated modes and wind-tunnel data 
points used in the prediction. The second volume presents the 
power spectral density plots made for all of the predictions. 



quite low when compared with the maneuver loads spectrumo Local 
loadings can be high, however, and can result in troublesome but 
not usually serious secondary structural fatigue. As a result, 
the prime concern with high intensity buffet is the dynamic 
response caused by buffeting loads and how it limits the opera- 
tional envelope of the total airplane system. 

The buffet problem, particularly transonic buffet, has 
received increased interest in government and airframe industry 
supported research programs. The Air Force Flight Dynamics 
Laboratory and the NASA Langley, Ames, and Flight Research Centers 
have supported many flight test, wind tunnel, and theoretical 
investigations aimed at developing the capability to better under- 
stand and predict buffet characteristics of new aircraft designs. 
Diverse opinions among the investigators has led to much confu- 
sion on the nature of the buffet problem and how characteristics 
may be predicted. All prediction methods currently in use are 
semi or purely empirical with which buffet onset can be predicted 
with fair accuracy o The problem of predicting buffet intensity 
is far more complex and requires consideration of many factors 
besides the aerodynamics of separated flow. 

A wind tunnel study was conducted at NASA Ames Research 
Center to determine the nature of fluctuating pressures during 
buffeting flow at subsonic and transonic speeds (ref. 1). The 
model used was a rigid 1/6 scale semispan model of the F-lllA. 

The wing had variable sweep and was heavily instrumented with 
high response pressure transducers. It was desired to correlate 
the model data with flight test results; however, this was no 
simple task. 

A supporting study was initiated at the Fort Worth Division 
of General Dynamics for NASA to collect and analyze flight test 
data for correlation and to develop a means for performing the 
correlations. Since the principle objective of the total program 
was to better understand the mechanics of high intensity buffet, 
the F-lllA was an ideal test bed because of the large amount of 
flight test data available from the loads program and the wide 
variety of Mach-altitude-wing sweep configurations available. 

Some results of analysis of this flight test data are given in 
references 2 and 3. The complete results of the flight test data 
analyzed are given in reference 4. These results were invaluable 
in the development of the prediction method to be given in this 
report. 
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This report dociiments the development and assessment of a 
method with which fluctuating pressure data obtained from rigid 
scaled wind-tunnel models can be used to predict flexible full- 
scale buffet response. The method requires unsteady aerod 3 mamic 
forces, natural airplane modes, and the measured pressure data 
as input. A gust response computer program is used to calculate 
buffet response due to the forcing function posed by the measured 
pressure data. This method is similar to that given in refer- 
ence 5 except that the latter does not include aerodynamic forces 
due to airplane rigid and elastic motions o By calculating both 
symmetric and antisymmetric solutions , upper and lower bounds on 
full-scale buffet response are formed. Final results are given 
in the form of ,upper and lower bounds on the power spectral densi- 
ties and the RMS values for angle of attack variation in maneuvers 
at several Mach-altitudes , Comparisons of predictions with flight 
test results are made and the effects of horizontal tail loads and 
static aeroelas ticity are shown. Discussions are also presented 
on the effects of primary wing torsion modes, chordwise and span- 
wise phase angles, and altitude. 

The report is broken up into two volumes. The first volume 
presents a description of the prediction method, how the method 
was evolved, and comparisons of predicted rms responses with 
flight test data for assessing the capability of the method. 
Included with the comparisons are complete descriptions of the 
flight test maneuver and data points, airplane conditions, and 
calculated modes and wind-tunnel data points used in the pre- 
diction. The second volume is composed primarily of the power 
spectral density plots made for all of the predictions. On these 
plots, the flight test power spectral density distributions are 
shown for a detailed comparison. Included with the plots are 
repeats of the flight test and prediction conditions for each 
case. 


SYMBOLS 


bj- reference length 

Fj. aerodynamic loading or moment due to 

the rth mode 

Fg aerodynamic loading or moment due to 

the buffeting pressure field 
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hj-(x,y) deflection of the rth mode at point x,y 

(Tjjj(x,y) mass concentrated at point x,y 


M 


free stream Mach number 


No 

q 


qr 


Qrs 


QrB 


a 

Ap(x,y, ft) ) 


6>i j ( ft) ) 

A 


ft) 


generalized mass in the modal coordinates 
for the airplane djmamic model 

characteristic frequency (frequency centroid 
of the power spectral density distribution) 

free stream dynamic pressure 

generalized coordinate 

generalized force in the rth mode due to 
aerod3mamics in the sth mode 

generalized force in the rth mode due to 
the buffeting pressure field 

wing angle of attack 

buffeting pressure amplitude at point x,y 
and frequency co 

phase angle between points i and j for the 
buffeting pressure field at frequency co 

wing sweep angle 

power spectra of the buffeting pressure 
field at point i and frequency ft) 

circular frequency 
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THE PREDICTION METHOD 


A method has been developed for predicting full scale buffet 
response characteristics by the use of wind tunnel buffet pressure 
measurements on rigid models and conventional gust response 
analysis methods. The wind tunnel pressure data is in the form 
of power spectral densities (PSDs) for each pressure transducer. 
Phase relationships (phase angles, ^ij) are available from cross 
power spectral densities (CPSDs) between pressure transducers. 

From the PSDs and phase angles, 0ij , complex pressure distribu- 
tions over the wing are constructed. The structural characteris- 
tics are obtained from a math model of the airplane with the 
proper mass distribution for the flight condition being analyzed. 
Aerodynamic stiffness and damping for the desired Mach-altitude 
condition are computed from oscillatory lifting surface theory 
based on linearized uniform potential flow. All of these items 
are combined into the d}mamic equations which are solved in the 
conventional manner of gust response analyses. 

The relationships between the above functions are shown in 
figure 1. Each block represents an individual computer program 
or computer run. The unsteady aerodynamic terms, [Qrs] and {Fi}^ 
are both computed with the same program. These terms are also 
expanded in frequency with a single programo Likewise, the 
buffet forcing function terms {Qig} and {^Fg) are calculated in 
the same mannero Each program used in the flow diagram will be 
described in the following subsections where major emphasis will 
be placed on the programs to convert wind tunnel data and to 
calculate the wing and tail pressures at the aerod 3 mamic panel 
points . 


Existing D 3 mamic Response Capability 

The dynamic response capability at the Fort Worth Division 
is composed of several computer programs, all of which together 
represent a state-of-the-art capability, typical of the major 
airframe companies. The basic programs are described in the 
following paragraphs. 

Computation of normal modes , - Normal modes of vibration are 
computed with a two-dimensional finite -element computer program 
with the Fort Worth Division designation of BY7 (ref. 6), The 
BY7 program computes a matrix of stiffness influence coefficients 
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for arbitrary vehicle configuration components such as wing 
surfaces, vertical and horizontal tail surfaces, and control 
surfaces. The procedure also has an anisotropic capability. 

When appropriate mass data is supplied, normal modes can be 
computed for cantilevered vehicle components or an entire 
vehicle with free-free boundary conditions. 

The main reasons for using the BY7 program rather than 
NASTRAN are economy and convenience „ Although NASTRAN has been 
set up at the Fort Worth Division and routinely used since 1970 
it has been found that for airfoil type surfaces, the cost in- 
creases by an order of magnitude with little or no gain in 
accuracy when compared with the BY7 programo Perhaps the 
strongest reason for not using NASTRAN is that a math model 
had been previously set up of the F-lllA/FB-lllA for BY7 
(refo 7). Tables 1 and 2 show the comparison between the pre- 
dicted frequencies obtained with the BY7 math model and those 
from ground vibration tests for Airplane #12. The total vehicle 
was represented with free-free boundary conditions. The close- 
ness of agreement clearly illustrates how well the model repre- 
sents the airplane. 

Mode interpolation and computation of generalized masses . - 
Mode interpolation for the unsteady aerodynamic programs is 
accomplished with a surface spline method which is applicable to 
entire airplane configurations. The method is based on the sur- 
face spline technique of Harder and Desmarais (refo 8) , The 
interpolation procedure also calculates the generalized mass 
terms, Mrs? with the following equation 

“rs " r"b 3 hi-(x,y) hg(x,y) <r^(x,y) d x dy 

^ ^ A 

where 

p - density 

bj- = reference length 

hj.(x,y) ,hg (x,y) = deflection of the r^^, s^^ modes 

at point x,y 

oj^(x,y) = mass distribution per unit area. 
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Computation of unsteady aerodynamics . - Several programs are 
available for computing unsteady aerod 3 mamics based on linear 
potential flow theory. These include kernel function methods 
as well as the doublet lattice and Mach box methods. 

The subsonic doublet lattice method (ref. 9) as it exists 
at the Fort Worth Division is designated as procedure TS5, This 
program has the capability of treating up to ten trapezoidal 
surfaces lying in the same plane. Thus, planar representations 
of wing -body -horizontal tail configurations can be achieved as 
has been demonstrated in previous work. As will be discussed 
later, this program is used to compute total airplane aero- 
d 5 niamics and make estimates of loads induced on the horizontal 
tail due to buffet loads on the wing. 

Supersonic aerod 5 niamics are calculated with a Mach box 
procedure based on the method of Pines and Dugundji (ref. 10). 
This procedure, designated as BQ8 at the Fort Worth Division, 
is limited to wing alone configurations. Hence, the supersonic 
response predictions in this report are limited to wing alone 
simulations o 

The aerodynamic programs are used to compute the generalized 
force matrices, [ Qrs ] • These terms are computed as 

“ 4obl o, 2 _/y'hr(x,y) [^-Apg(x,y, w )] dxdy 
A„ 

where 


p = density 

bj- = reference length 

0) = exciting frequency 

hj-(x,y) = deflection of the r^h mode at point x,y 

~ pressure distribution amplitude induced 
at point x,y by the sth mode oscillating 
at frequency, co 

= wing area 
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other terms that are computed are the Fr terms which represent 
aerodynamic loads, moments, and stresses induced by each mode for 
response calculations* These terms will be discussed in the 
section on d}nnLamic response. 

Computation of [Qrsl at additional frequencies . - 

The dependence of the [Qrs] and -{Fr} terms on frequency is 
implicit and very complex. The variation with frequency is 
generally quite smooth, however, and is easily amendable to 
interpolation or approximation schero.es with a high degree of 
confidence. 

An approximation scheme has been developed at the Fort Worth 
Division which has proven to be very dependable and efficient 
(ref. 11). The method uses a parabolic fit to condition the 
terms near zero frequency. The conditioned terms are then fitted 
with Tschebychev polynomials of the first kind for up to 10 
fitting points. By the use of an exponential transformation, 
the approximation becomes theoretically valid for the frequency 
range 0 < oj < oo , In practice, the upper and lower values input 
are used as the real limits. 

Shown in figures 2 and 3 are some examples of application 
of the approximation scheme. The Qn term shown in figure 2 
is the real part generalized force resulting from a rectangular 
wing of aspect ratio 2 oscillating in a unit translation mode 
for M=1.2. The circles are the fitting points and the squares 
are values computed at other points » As can be seen, the 
approximation is quite good at off-fitting points. In figure 3, 
the real part of Ql2 is shown where the second mode is wing pitch 
about its leading edge. With the approximation fitted at the 
circle s 3 mibols, the agreement with the square symbols is again 
quite good. 

This method was used to expand the [Qrsl and '{Fr} terms as 
calculated with the unsteady aerod 3 mamic programs. The calcu- 
lated values were obtained at six fitting points and expanded to 
match the full scale frequencies of the wind tunnel model data. 

The approximation method is available in the doublet lattice 
procedure, TS5, and is usable as a separate capability if desired. 
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Computation of dynamic response .- The d 5 rnamic response 
procedure, designated as SB2 at the Fort Worth Division, solves 
the equations of motion for such linear systems as a total air- 
plane subjected to an arbitrary harmonic forcing function 
(ref. 12) o The airplane must be balanced aefod 3 mamically, thus, 
rigid body motions are obtained as well as elastic response 
motions. The equations of motion are expressed in terms of 
generalized coordinates. In matrix notation, this becomes 


^1 

! 

All Ai 2 . . . AiHq j 

-1 

[ 

^2 


! ; 

A 21 A 22 • • • A2Nq 

i 

Q2f 

< 

• 1 

• j 

1 -- 

• « • 1 

• • • 

< 

j 

i 

i 

• f 

• 

! i 

i 

1 ! 

1 1 

0 » • 

AlSIqi Ajqq^ . . . Aj^q^q 

1 1 

0 i 

L 


or 

f— r -| “1 r 'I 

< ^r > ~ ^1 ^^f ^ 

L J L -J 

where 

Arr ^ (1 + i Sr) - i2rr (■'^) ] ^^rr + Qrr. r=s 

Ars = Mj-s + Qrs, 

(0 ~ exciting frequency, rad /sec 

= undamped natural frequency of the r^^ 
normal mode, rad /sec 

gj- = structural damping coefficient for the r^^ 
mode (usually assumed as gj- = 0.03) 

7j- = ratio of viscous damping to critical damping 

Nq - total number of generalized coordinates or 
degrees of freedom, l<Nq<20 
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M]-s “ generalized mass as described previously 

Qj-g = generalized work of the unsteady aerodynamic loads 
on the oscillating airplane as described previously 

Qj~f = generalized work of the forcing function on the 
oscillating airplane 

i = 

The Q 2 :f terms are computed in the same manner as the Qj-g terms 
where Apf(x,y) of the forcing function is used in the integral 
rather than Aps(x,y). 

Once the equation of motion is solved for each frequency of 
the forcing function, the ‘{qr}* ^.re obtained with which deflec- 
tions, velocities, accelerations, loads, moments, or s'tresses 
can be computed. 

Normally, in procedure SB2, the response characteristics 
are computed as transfer functions in terms of response per unit 
input as a function of frequency o Total response is then calcu- 
lated with an input spectrum of the forcing function (a Dryden 
or Von Karman gust spectrum, etc.) to yield an output spectrum. 
The transfer function, HL(<i >)5 for the Lth point and exciting 
frequency, o) , is expressed as 

Nq _ 

Hl( (^) = [ (f " f’rL + ^rh( ^ ) J 

r=l 


+ (i w FfL + FfL< oj ) 

where 

Nq = number of degrees of freedom or 
modes (1< Nq^ 20) 

(i o) )^FrL ~ inertia component of the response in 
each mode where p is a function of 
the type of response item 
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Ftl( ^ ) “ load, moment, or stress component of 
the response in the rth mode 

(i<a)PFfL “ inertia component of the response due 
to the forcing function 

FfL( m ) = load, moment, or stress component of 

the response due to the forcing function 

The above terms will be discussed in more detail. 

The first term, (ia))^FrLj which is the contribution of 
inertia to the response point, will vary in form depending 

on the type of response desired. If acceleration response is 
desired at the point, then p - 2 and Fj-^ = h^^j^ or 

(i Fj-l = (i c»)^ hj-L 

where h^L Is the deflection of the r^h mode at point L. 

If velocity or displacement is desired then p = 1 or 0, 
respectively. For loads computation, p = 2, and Fj^l defined 
for shear as 


r(ia.)VLl 

^ SHEAR 



m(x,y)hr(x,y) dxdy 


where 


m(x,y) = mass distribution 
hj-(x,y) = rl-^ mode shape 

Al = area for calculating the total shear 

The area, Al, is defined by the planform and the station and 
axis at which the shear is desired. For bending moment, FrL 
becomes 

[(ift))^FrL] ni(x,y)hj.(x,y) dxdy 

B.M. 
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where (Y-Yref)'' is the moment arm perpendicular to the bending 
moment axis. Likewise, torsion is defined as 

r(i(u)^FrL] = - ff (x,y) dxdy 

TORS. 

where again (X-Xj^p) is the moment arm perpendicular to the 
torsion axis. 

The second term, FrL(^)j is simply the total aerodynamic 
load, moment, or stress induced at the Ltk point due to the 
motion of the r^k mode. For the computation of displacements, 
velocities and accelerations, FrhC ^ ) is set to zero. It is 
expressed as 


FrL^ ^ ) = // Apr(>!^,y) 

Al 

where 

Apr( x^y) ~ unsteady aerodynamic pressure distribution 
induced by the r^i^ mode of unit amplitude 

^ > shear 

= (Y-Yref) ' > bending moment 
= (X-Xref)* j torsion 


The terms for f^(x,y) are the same as those used in the compu- 
tation of Fj-l. 

The third term, (iw)^FfLj i® inertia component which 

is due to the forcing function. For gust response, this term is 
zero, hence it is also zero for buffet response. Cases for which 
the term is non-zero are response due to forced control oscilla- 
tion, forced fuel sloshing, or internal inertial shakers for 
flutter testing. The expression for (iw)PFfL identical to 
that for (i o) )PFj-l- 

The fourth term, FfL( ) > is the aerod3mamic load, moment 
or stress induced at the Lth point by the forcing function. 

This term is also zero if displacements, velocities, or accelera- 
tions are being computed o The expression is 
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FfL<") = // Apf(x,y) fL(x,y) dxdy 

Al 

where Apf(x,y) is the loading distribution of the forcing 
function and f^, is as defined previously. 

For the present task of calculating buffet response, the 
input forcing function is actually "shaped" as will be shown in 
the next section. Thus, the transfer function as computed from 
the above quantities will not be a true transfer function but 
will instead be the square root of the power spectral density. 
This is seen from the expression for calculating the output 
spectrum, 02(<wi)j resulting from the effect of a transfer 
function, H2l(oJi) on an input spectrum, 

02(a>i) = [H2i(wi)j^ 01 (wi) 

In the case of buffet response, the 0i is set to a constant 
value of unity, i.Oo, 

where coi and are the first and last frequencies in the 
spectrum being analyzed. 

Results from SB2 are in the form of a PSD, total RMS, and 
No of each response item. Other data are also output, however, 
they are not of major interest at this points The RMS values 
are calculated over the frequency range of the prediction, hence 
comparisons with flight test data must be made over the same 
frequency range o The No values are also calculated over the 
same range o This term represents the total number of times that 
the function being analyzed crosses the mean axis with a positive 
slope. Hence, No is a characteristic frequency. Where RMS is 
calculated as 



RMS 



No is calculated as 


No = r (Aco)± 0 ( a>i) 

rS 

In addition to the PSD, the integrated quantities, RMS and No, 
provide a convenient means for comparison with flight test. 

The next section will describe the process by which the 
wind tunnel pressure PSD and CPSD data are £ombined with the 
structural modes to produce the {^Qrf}* and *(Fr^ terms used in 
the d 3 mamic response programo In order to distinguish between 
the more conventional forcing function terms and those due to 
buffeting pressures, these terms are designated as and {Fg}, 


The Buffet Forcing Function 

Referring to figure 1, the process of converting wind tunnel 
data PSDs and phase angles into generalized forces is a three 
step process. First the wind tunnel data is transformed into a 
complex pressure distribution over the wing for each frequency 
at which a PSD estimate is made. These pressures are then inter- 
polated at the load points on the doublet lattice panels. The 
interpolated values are used in conjunction with the doublet 
lattice unsteady aerod 5 rnamic matrices to calculate estimated 
pressures induced on the horizontal tail by downwash produced 
by the wing buffet pressures o Finally, the wing and tail pres- 
sures are used in the doublet lattice program to calculate the 
generalised forcing function terms, j the response 

terms, <^Fg}^ , for shear, bending moment and torsion on the wing 
and tai 1 . 

The following paragraphs will describe (1) the means by 
which static aeroelastic effects are accounted for, (2) conver- 
sion of the wind tunnel data, (3) calculation of horizontal tail 
buffet pressures, and (4) calculation of the generalized forcing 
function and response terms. 

Calculation of static aeroelastic effects .- The wind tunnel 
buffet pressure data is obtained from a rigid model and hence 
does not represent the flexible airplane in flight. In order to 
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approximately account for the static aeroelastic effects, the 
rigid angle of attack values, orig, used to obtain the wind 
tunnel data are modified with theoretical results taken from 
aeroelastic analyses of the F-lllA. The appropriate flexible 
to rigid lift coefficient ratio is used to calculate the angle 
of attack, apLEX? required for the flexible airplane to produce 
the equivalent lift that would be produced by the rigid airplane 
at aRiG. 

The flexible to rigid ratios used in this study are summa- 
rized in figure 4 for the three wing sweeps considered. The 
values for apLEX ^Lre calculated as 

^ ^ RIG 

pL flex X 

RIG M, h, A 

for the desired Mach-altitude-wing sweep condition. Thus, the 
buffet response calculated for a given wind tunnel model angle 
of attack, a rig? is used as a prediction for an airplane angle 
of attack, apLRx- 


This technique is not exact in that it assumes that the 
aeroelastic twist has no effect on changing the buffet pressure 
distributions, that is, the buffeting forces are directly related 
to total normal forces on the wing. It is suspected that aero- 
elastic twist tends to "soften” the buffeting pressures which 
makes the above assumption conservative o If the suspicion is 
correct, then even the predictions which are shifted in a due 
to gross aeroelastic effects should overpredict flight test 
results. Since the results which will be presented later in 
this report are overpredicted, it is felt that this reasoning 
may be valid, however, further research is needed in this area. 

Conversion of wind tunnel data to pressure distributions . - 
The wind tunnel data which is used to generate buffet forcing 
functions in the prediction method is received on magnetic tape 
as a collection of PSDs , CPSDs , phase angles, coherency functions 
and convection velocities. The data is available for multiple 
sets of three pressure transducers which are referred to as x, y 
and Zo A complete coverage of the transducers on the wing is 
available for each Mach and angle of attack condition. The 
means by which this data is processed into complex pressure 


distributions over the wing is embodied in a computer program, 
the methodology of which will be described below. 

The only data that are used on the magnetic tape are the 
PSDs, 0m(t«i) j phase angles, 6mn(c^i)- These quantities 

are combined to obtain complex pressure distributions in the 
manner illustrated in the following example for a chordwise 
distribution. Consider the two dimensional distribution in 
figure 5a where six transducers are located along the chord. 
Because the data processing used at NASA ARC on the wind tunnel 
data requires that the data be composed of sets of three items, 
three sets are used to describe the 0m(^^i) and 0mn(wi) for the 
example as shown « The first set translates to 

= (01 . 02 > 03 ) 

(exy. ^xz » ®yz )l = (®12> ®13» ®23) 
the second set to 

(0x j 0y j 0z >2 " (^3 > 04 » ^^5 ) 

(®xy’ ® xz ’ ®yz '>2 “ ®35> ®4S) 

and the third set to 

(0x 5 0y j 0z )3 = (05 > 06 » 07 ) 

(®xy’ ®xz> ^yz )3 (®56> ®57> 

where each term is a function of frequency. The 07 may be a 
repeat of 05 or some other quantity which is not needed. 

The main problem of translating the above data into a 
complex pressure function is determining the reference point 
from which the phase angle is measured „ Currently, a point is 
chosen as the "leading edge" point and its phase angle is set 

at zero. Points which are forward of that location are nulled 

to zero in both phase and amplitude. In this manner the effect 
of shock induced separation can be simulated. The determination 
of the forward reference point is based on a scan of the PSDs 
and the coherency functions along the chord by the program user. 
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The coherency functions are usually low or show considerable 
scatter for points forward of the shock or the point of separation. 
Aft of the shock or separation point, the opposite is true which 
indicates good correlation particularly for adjacent transducers. 

For the example, let the reference be the second transducer. 
Thus the amplitudes become 


01 = 

= 0.0 


= 0.0 




02 = 

= 0 

02 

= 0.0 




03 = 

= 03 

®3 

= 023 




04 = 

= 04 

04 

= 023 


034 


05 = 

= 05 

05 

= 023 

+ 

034 + 

045 

06 = 

= 06 

06 

= 023 


034 + 

045 +056 


ignored Gy 


and the complex pressure distribution is obtained as follows 

= 0.0 

q\/(A")i 

[ 02 (" i)J ^ 

[ 03 ("i)] ^ EXPrie3(6)i) 

[06(<"i)] ^ EXp[i06(<»i)] 


Ap2(g)j[) 

qy(Ao>)i 

Ap3(<ai) 

q>/(At«)i 

q7(Ao>)i 
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where the dependence on ct)± is indicated for clarity. The (Ao>)i 
quantity is the finite bandwidth used in the NASA ARC spectral 
analysis. Its value changes as o)i in order to provide a wide 
frequency coverage with a minimum of estimates. 

Three important characteristics of the above scheme should 
be noted. First, the only phase angles used are those between 
adjacent transducers, i.e,, only 0xy or 0yz are used. These are 
used primarily because the coherency between adjacent points is 
better which leads to a higher confidence in the cross correla- 
tions. Second, the items to be nulled are still used but set to 
zero and those to be ignored are discarded. Thus care must be 
exercised in the selection of points to i nsure that all points 
are used properly. And third, the \/ (A&)) i quantity is retained 
in the pressure distributions. This is done so that the result- 
ing response transfer functions will be of the form /PSD which 
eliminates points of confusion that would be caused by a varying 
(A o>) i . 

The only difference between the two-dimensional case de- 
scribed above and a three dimensional case for an entire wing is 
that several chords must be treated and the spanwise phase angle 
must be accounted for. In each chord a reference point is selec- 
ted and a set of transducers is assigned. The phase angle con- 
struction then proceeds for each chord in the manner described 
above. The spanwise data is received in the same form of 0 and 0 
for the running loads o These data are processed in the same 
manner as the chordwise data into the form of complex running 
load distributions. The reference point in the spanwise data is 
always taken to be the inboard chord station which has a zero 
phase angle. The spanwise data are used to determine the proper 
absolute phase angles for the chordwise data. 

The spanwise data are physically calculated as integrals of 
the chordwise data at each span station. Thus, the phase angle 
of the integral of the processed chordwise pressure must be equal 
to that of the processed running loads data at its respective 
span station. Rather than calculate the phase angle change 
necessary to satisfy this condition at each span station, the 
complex chordwise pressures are processed directly as 

Apnm(<^i) = ^Pnm(^^i) ^£m 

_q>/( qy( Aft))i jC^nii j/iPm I 
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where 


= complex pressure value at the 

chordwise point on the span station 

with the desired absolute phase angle 


= complex pressure value with phase angle 
referenced to zero at the leading edge 
independent of span station 

= complex running load at the m*-^ 
span station 

= complex conjugate of the chordwise 
integral of Apmn(^i) 

q^( Aco)i 

The final result from the wind tunnel data conversion 
program is a magnetic tape containing the following data; 

ft>^br 

(1) The reduced frequencies, for each 

spectral estimate, k]_, ... k^. 

(2) Complex pressure distribution, 

qsJ(Aco)^ 


, ’ ^ . Ap(x,y,ft>T^) 

(K + 1) Complex pressure di.stri.butr on, — — 

q/(Ao))j^ 

For K total frequencies. This tape is valid for a single Mach- 
alpha condition assuming that density and/or Reynold’s number 
effects are of second order magnitude. 

Calculation of the horizontal tail buffet pressures . - 
Horizontal tail buffet pressures were not measured in the wind 
tunnel test program (ref. 1). However, in the prediction of 
airplane buffet response, these data were needed for the forcing 
function. Thus, a method was developed with which the tail 
pressures could be estimated with linear theory aerodjmamic 


Apnm(^i) 

qy< Aa>)i 


0 . 


^Pnm(^i) 

q^( A aOi 


'£m 


Ap 


19 



matrices. It was assiomed that the wake in the vicinity of the 
horizontal tail due to buffeting pressures on the wing could be 
predicted with doublet lattice unsteady aerodynamic influence 
coefficients from a known wing load distribution. Such an as sump 
tion is not unreasonable for distances of several panel chords 
downstream. 


In matrix form, the aerodynamic problem is 

[ Aww Awt"| r pwl = f Ww"' 

Atw AttJ \ PtJ LWt/ 

where the aerodynamic influence coefficient submatrices are 


Aww = wing on wing 
Awt = tail on wing 
Atw == wing on tail 
Att ~ tail on tail 


and 


{pt} , {Wt} 


pressure and downwash on the wing 
pressure and downwash on the tail 


Since the wing is buffeting, pw is known from the wind-tunnel 
data as pwB- T^ie horizontal tail is usually at a negative angle 
of attack during a high- a maneuver; hence, it is assumed to be 
in a buffet free condition. Also, the tail upper surface is 
effectively the only portion feeling the buffeting wake from the 
wing upper surface; thus, 

%[Atw] {pwb} + [Att] {PtB} = {wt} = 0 
which yields 

{Ptfi} = [Att]' [Atw] {pwb} 


With this form, all that is necessary to obtain ptg from pwB is 
to calculate the doublet lattice unsteady aerod 3 mamic matrices 
(ref. 9) at the same frequencies for which wind tunnel data PSD 
estimates are made. In actual practice, the matrices are not 
calculated directly at all frequencies but are interpolated from 
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a smaller set since about forty frequencies are used. Otherwise 
the cost would be prohibitive. 

Calculation of the buffet forcing function and response 
terms . - The doublet lattice program (ref. 9) is used to combine 
the buffet pressures with structural modes to obtain the general- 
ized forcing function terms. It is also used to determine the 
shear, bending moment and torsion due to the buffet pressures for 
response calculations. The aerod 3 naamic matrices generated by the 
doublet lattice program for the wing- tail unsteady aerod 3 mamic 
pressures are also used to predict horizontal tail pressures due 
to the downwa'sh produced by the wing buffet pressures in subsonic 
flow. Thus, the first step in obtaining these data items is to 
interpolate the wing buffet pressures from the transducer loca- 
tions to the doublet lattice panel load points. 

The interpolation of the buffet pressures is performed with 
a surface spline technique similar to that used to interpolate 
the structural mode shapes. The surface spline simulates the 
deflected shape of an infinite plate pinned at the points which 
are being interpolated fromo The governing equation is 

N 

bn(rnj)^ ln|( Hj)^| 

n=l 

For n=l, N where wj is the function value at point (xj,yj) and 

+ (yn-yj)^ 

The points (xn,yn) the fitting points which in this case are 

the transducer locations. The coefficients [ai, a2, ^ 3 , bi,...bNJ 
are determined by equating the wj to the known pressures at points 
(xj,yj) = (xn,yn) n=l, N. Three other equations are also 

satisfied, 

N N N 

XI bn = xn bn = X^ 7n bn = 0 

n=l n=l n=l 

to give the necessary N+3 total equations to solve for the N+3 
unknown coefficients. With the known coefficients, the pressures 
at the doublet lattice points are calculated. 
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An example of the surface spline fit is shown in figure 5b . 

The condition is M=0.80, aw=10.18O ( abody=9. 18°) and fmodel=48 Hz. 
The two span stations shown are ^=0.600 and 'i?=0.77* The plotted 
S 3 mibols are the processed values obtained from the wind tunnel 
data conversion program. The surface spline is the smooth curve 
which shows a well behaved fit of the data. This comparison is 
typical of other conditions and serves to illustrate the power of 
the surface spline to smooth the incredibly ragged data which is 
being used. 

The wing pressure values obtained from interpolation, PwBj 
are the values used to calculate the tail pressures, ptB, 
discussed in the previous section. The pwB 3,nd PtB ^.re then 
input to the^doublet lattice program for final calculation of 
the QrB and Fb terms to be used in the response procedure, SB2, 


Buffet Response Prediction 

With the input data completed for the d 5 mamic response 
program, the final step in figure 1 is to calculate the buffet 
response in terms of accelerations, loads, and moments. Pre- 
dictions are made for both pure S3mimetric and antis5mimetric 
airplane motions. Since flight test results (refs, 2, 3, 4) 
indicate that the airplane response is usually as 3 mimetric, even 
in a "pure" S3nnmetric maneuver, both the symmetric and anti- 
symmetric responses are combined to produce an upper and lower 
bounds on buffet characteristics. These bounds are given as a 
function of angle of attack at a particular Mach-altitude-sweep- 
gross weight configuration. Since airplane buffet is subject to 
many variables other than those already accounted for (such as 
pilot characteristics and atmospheric conditions) the concept of 
a predicted bounds seems to be a very attractive means for treat- 
ing the high intensity buffet problem. 

The upper and lower bounds spectra, 0u a,nd 0 l> ^.re based on 
the following assumptions: 

1, The maximum response possible is obtained when 
both S3nnmetric response spectra, 0s, and anti- 
symmetric spectra, 0 a, are in phase and 1007o 
active at all frequencies: 

0U =[0s^ + 0A^] ^ 
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o The minimum response possible is obtained when 

(1) 0s and 0A are 100% active and 180° out of phase, 

(2) 0s is active only, or (3) 0 a is active only: 



The 0u and 0 l can also be used to obtain upper and lower bounds 
on the RMS values of response. Examples of both will be presented 
in the application of the method. 

The above equations for 0u and 0 l represent two extremes 
between which all flight test results should fall. Since 100% 
excitation is not possible, it would be expected that excursions 
outside of the bounds would be more frequent for the lower than 
the upper. Hence, the norm of the flight test data for all cases 
should be weighted toward the lower bounds. Comparisons between 
predictions and flight results confirm this hypothesis. 


EVOLVMENT OF THE PREDICTION METHOD 


Throughout the development of the prediction method, 
extensive comparisons with flight test data were made in order 
to determine the importance of various phenomenon. Results of 
these comparisons were used to guide the development of the final 
form of the prediction method which was summarized in the previous 
section. This section presents a summary of the more important 
effects as observed during the method development. 


Effects of Static Aeroelas ticity 

Shown in figure 6 are the predicted upper and lower bounds 
RMS values for the wing bending moment and the corresponding 
characteristic frequencies, Nq- Also shown are the flight test 
results at the same configuration (given in the figure) . The Nq 
is the frequency centroid of the PSD curve; hence, its use in 
conjunction with the RMS provides a much better means for compar 
ing integrated quantities. It is. also interesting to note the 
variation of No with ot in the results to be presented. 
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Figure 6 illustrates the improvement of agreement with flight 
test results if the predicted a values are shifted according to 
the flexible to rigid lift coefficient ratio as determined from 
figure 4. The plots on the left are shown with the predictions 
plotted at the rigid values of a from the wind tunnel data. The 
right hand plots use the flexibilized value of o«. In both cases 
the flight test data are plotted at their actual a. values. 

All of the predictions to be given in the following paragraphs 
will be shown as functions of the flexible a. values. 


Asymmetric Responses 

The examination of flight buffet data has clearly shown that 
in general the airplane responds in an asymmetric manner under 
high intensity buffet (refs. 2 and 3). The presence of major 
response in both symmetric and antisymmetric modes led to the 
formulation of the upper and lower bounds concept as the only 
reasonable means by which airplane buffet response could be pre- 
dicted, Thus, such uncontrollable items as pilot characteristics, 
atmospheric turbulence, airplane mass and geometric asymmetries, 
and control system inputs could be covered to a large extent. 

The upper bound was taken to be the sum of both predicted 
symmetric and antisymmetric responses acting in phase at all 
frequencies. Although this may seem conservative, it will be 
shox^ that this point can sometimes be exceeded. The lower bound 
was taken as either the symmetric or antisymmetric responses or 
the differences between the symmetric and antisymmetric responses, 
whichever was the lowest value at that corresponding frequency. 

An example of an upper and lower bounds PSD plot is shown in 
figure 7 for the wing tip accelerometer. Comparison is made with 
both the right and left wing tip results from flight test. It will 
be noted that even though there is a wide separation between the 
bounds, there are several points at which the flight test data 
confirms the difference. At about 12, 21, and 27 Hz, the right 
and left wing tips indicate such a superpositioning of symmetric 
and antisymmetric responses. Also, exceedence of the upper bounds 
from about 15 to 21 Hz was felt to be due to some non-linearities 
not covered in the prediction method. 

Another example of the upper and lower bounds is shown in 
figure 8 to illustrate variation with a. This prediction is for 
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the wing bending moment (same condition as in figure 7) where 
comparison is made with one wing only. These plots are the PSD*s 
for the RMS -No results given in figure 6. Comparison of figures 6 
and 8 shows how the change in spectral shape toward the first wing 
bending modes (s 5 mi. = 4.8 Hz, anti. = 7.4 Hz) with increasing 
buffet intensity is reflected by a steady decrease in Nq. The 
flight test data are well bounded by the predictions. 


Effect of Horizontal Tail Loads 

During the prediction method development the effect of adding 
horizontal tail unsteady aerodynamics and buffet loads was studied. 
Shortcomings of the initial prediction method which used wing aero- 
dynamics only, led to this investigation. However, measured buffet 
pressure data were not available for the horizontal tail since the 
model was not instrumented as such. Thus, the semi -empirical 
scheme which used the doublet lattice aerod 3 raamic matrices as 
described in the previous section was developed. 

The results shown in figure 9 illustrate the effect of vary- 
ing horizontal tail loads on the wing shear. The wing alone 
results are obtained with wing buffet pressures and wing unsteady 
aerodynamics. For this simulation, the lower bounds seem to be 
more representative of the flight data and the No plot does not 
agree at all. The PSD*s for this case verified the disagreement. 
The second solution with total airplane unsteady aerodynamics and 
full horizontal tail buffet loads showed a significant improvement, 
particularly for the Nq comparison as was reflected by the PSD's. 

It was apparent, however, from examination of the horizontal tail 
pivot loads, from both flight test and prediction that the tail 
buffet pressures were too high. Further thought on the matter led 
to the conclusion that the estimated tail buffet pressures should 
be divided by two. This conclusion was based on the fact that at 
high angle of attack, the buffet wake which leaves primarily from 
the wing upper surface has very little effect on the tail lower 
surface. In addition, the tail is usually at negative incidence 
relative to the wing. The final solution with 1/2 horizontal tail 
buffet loads as shown in figure 9 verifies the reasoning in that 
the RMS loads are well bounded and the No results show excellent 
agreement. (This solution has another variation included where 
the wing first torsion frequencies of the airplane and wind tunnel 
model were matched which reduced the frequency limit of the pre- 
diction from 38 to 31 Hz. The frequency matching contribution to 
the improved agreement was insignificant compared to that due to 
using 1/2 tail buffet loads.) 
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Figure 10 shows the effect of the same variation of horizon- 
tal tail buffet loads on the wing tip accelerometer predictions. 
For this response item, the final method does not appear to be 
any more accurate than the wing alone method. Since the wing tip 
accelerometer is sensitive to all wing modes and the wing root 
shear is more affected by the lower wing bending modes this is 
not surprising. Total airplane motion which is affected by the 
horizontal tail has a greater influence in the lower wing mode 
responses, hence it would be expected to significantly affect the 
wing root loads. Likewise, it would not have as great of an 
influence on the higher wing modes which are equally important 
for wing tip motion. 

Figure 11 shows the results for the C.G. vertical acceler- 
ometer. Although the wing alone results show excellent agreement 
for a = 6.9° and 11.7®, they do not have the decrease in response 
at a = 14.1° as indicated by the flight test data. Also, since 
it is felt that it is impossible for the airplane to response at 
100% as the wing alone results show, the final method seems to be 
more realistic. 

From this study it was concluded that the horizontal tail was 
important in the prediction of airplane buffet characteristics 
under conditions well beyond buffet onset. It also appeared that 
the concept for estimating the tail buffet loads was correct; 
however, even with the 1/2 factor, the predicted tail loads were 
high as compared with flight test data. 


Effect of Buffet Pressure Spanwise and Chordwise Phasing 

The effect of including various combinations of spanwise and 
chordwise phase angle data in the buffet pressures used to form 
the forcing function were studied. The variations included use of 
the span loadings with and without spanwise phase angles as well 
as chordwise distributions with and without spanwise phase angles. 
The spanwise loadings were assumed to have a constant chordwise 
distribution, hence, generalized forces resulting from torsion 
modes or those with significant twist were not accurate for the 
variations which made use of these loadings. 

Shown in figure 12 are the PSD results for the wing tip 
accelerometer for the case discussed in the previous paragraphs. 
Predicted responses are shown for symmetric conditions only. 
Figures 12a, 12b and 12c correspond to the four variations made 
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for a — 6 . 6 ^, 11.1° and 14.4°, respectively. The buffet loading 
variations noted in the figures are defined as: 

. No phasing - Spanwise loadings without spanwise 
phase angle data as obtained from the processed 
wind tunnel buffet data (section normal force 
data discussed in ref. 1); 

. Spanwise Phasing - Spanwise loadings with spanwise 
phase angle data; 

. Chordwise Phasing - Chordwise distributions with 
chordwise phase angle data and leading edge 
pressures assumed to be in phase at all span 
stations; and 

. Spanwise and Chordwise Phasing - Chordwise 

distributions with chordwise phase angle data and 
spanwise phasing as determined from the spanwise 
loading data. 

The techniques with which the buffet pressure data are used to 
form the forcing function have been discussed under the paragraph 
''Conversion of Wind Tunnel Data to Pressure Distributions." 

Comparisons with flight test data made in figures 12 are 
given only for relative purposes since the predictions are made 
only for S 3 mimetric responses. It can be seen in all cases, 
however, that the inclusion of the elastic degrees of freedom 
(9 DOF) in the predictions result in consistently improved solu- 
tions as compared to those obtained with rigid body translation 
(1 DOF). Although individual peaks are not directly comparable, 
the order of magnitude of the predicted flexible PSDs agrees very 
well with that of the flight test data. 

Comparisons of the four variations in figures 12 for the 
1 DOF response curves show the direct effect of phase angle 
variation on the total fluctuating normal force. Since this 
degree of freedom is vertical translation, the total normal force 
PSD is directly proportional to the wing tip accelerometer PSD. 

In all cases the largest total normal force results from the use 
of no phasing as might be expected. The inclusion of either span- 
wise or chordwise phasing or both reduces the total normal force 
in varying degrees. The most notable change is seen in figure 12c 
where inclusion of chordwise phasing shows little change in the 
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1 DOF curve but inclusion of spanwise phasing results in a drastic 
change. In this case, the flow has become separated at the lead- 
ing edge and the buffet pressures apparently have such large span- 
wise phase angles that the total normal force becomes nearly zero 
at 15, 25 and 32 Hz. This canceling effect is noticeable in 
figure 12b (a= 11.1°) to a lesser degree at different frequencies 
and even less in figure 12a (a= 6.6*^). As a result it could be 
concluded that spanwise phase angle is more important for vortex 
or leading edge separated flows than for conditions at lower 
angles of attack where the buffeting forces are due primarily to 
oscillating shocks or shock induced separation. 

Effects of the four variations on the flexible airplane 
responses vary with the value of a and the type of mode shape. 

The first wing bending mode (4.794 Hz) shows the greatest response 
for no phasing and least response for both chordwise and spanwise 
phasing f or a = 6,6° anda= 11.1°. Fora= 14.4°, however, the 
response for spanwise and chordwise phasing is almost as large as 
no phasing. The first wing torsion mode (24.064 Hz) is most 
sensitive to chordwise phasing which is well illustrated in 
figure 12a (a = 6,6°). At the higher angles, the effect of phas- 
ing is not too dramatic. This is felt to be more due to the 
nature of the response item since, as will be shown, wing torsion 
response is sensitive to chordwise phasing at all angles _of attack 
Finally, a combined wing second bending mode with some twist 
(37.573 Hz) always shows maximum response for no phasing and 
minimum for spanwise and chordwise phasing with a range of almost 
two orders of magnitude f or a = 14.4°. This mode shows greater 
sensitivity to spanwise phasing at a- 11.1° and 14.4°, and about 
the same sensitivity to either chordwise or spanwise at a= 6.6°. 

Figure 13 shows the effect of phasing on the total RMS and Nq 
results for various accelerometers and wing loads and moments. 
Almost all items have the highest response for no phasing at all 
angles of attack with exception of wing bending moment ata= 14.4° 
The minimum response varies with angle of attack where chordwise 
and spanwise phasing is minimum ata= 11.1° and spanwise phasing 
is minimum ata= 14.4°. 


The Influence of Wing Torsional Motion 

The strong coupling between torsional motion of the airfoil 
and shock oscillation is one of the important characteristics of 
transonic Mach buffet. As the airfoil oscillates in low frequency 
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torsion, the shock on the upper surface moves forward with nose up 
displacement and aft with nose down. Thus, the lift increment on 
the airfoil due to shock movement is negative for a positive incre- 
ment in angle of attack, Aa. For low values of mean angle of 
attack, oe , the total incremental lift due to A a is positive as 
a result of increased lift forward of the shock. As a increases, 
the mean shock position moves forward on the upper surface in an 
orderly fashion until it nears the crest of the airfoil relative 
to the frees tream direction. As it nears the crest, it moves 
larger distances for a constant Aa, and, at the crest, it becomes 
very unstable and tends to move forward of the crest for a very 
small Aa. At this point, the total incremental lift due to the 
small A a is negative in two-dimensional flow as a result of shock 
movement. The additional lift upstream of the shock due to the 
small Aa is negligible; hence a "glich" results in the lift-vs.- 
a curve . 

In a dynamic sense, the above mechanism creates a very strong 
forcing function, which results in "significant buffet response at 
the first wing-torsion-mode frequency. This "torsion mode peak" 
should occur at a « 10*^ for Mach = 0.80 and A = 26° on the rigid 
F-lllA, as indicated by the rigid-body responses in the predicted 
PSDs. Approximate accounting for static aeroelastic effects 
increases the value to a « 10.9° at an altitude of 6035M and 
to a « 11.5° at 1490 m. The increased twist of the wing due to 
static aeroelas ticity , however, will tend to soften the mechanism 
by preventing the shock from crossing the crest at all span sta- 
tions at the same value of a . This effect cannot be accounted 
for in the linear static aeroelastic correction but is undoubtedly 
significant. The significance has long been observed as a 
flattening or complete disappearance of the transonic "glich" 
on static aeroelastic models. 

Shown in figure 14 are comparisons of predicted rigid and 
elastic responses with flight test results for several variations 
of the theoretical model. In figure 14a results are shown for a 
wing-alone prediction with directly scaled wind tunnel model data. 
These results show how the "torsion mode peak" in the model data 
at a ~ 10° (rigid) did not scale to correspond with the airplane 
first torsion mode frequency. In reality, this peak was scaled 
to occur at about 30 Hz as can be observed in figure 14a rather 
than at about 25 Hz which is the airplane first torsion frequency. 

Figure 14b shows how prediction compares with flight test 
results for a theoretical model which uses total airplane aero- 
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d 3 mamics and matched model and airplane first torsion frequencies. 
The total airplane aerod 3 niamics included 100% buffet forces on the 
horizontal tail. The frequencies were matched by changing the 
frequency scale factor such that the model first wing torsion mode 
frequency became 25 Hz. Comparison of figures 14a and 14b shows 
that although frequency matching resulted in improving the pre- 
diction in the vicinity of 25 Hz, the large peak observed in the 
flight test data was still not predicted. 

Figure 14c shows results of the final perturbation where 1/2 
of the horizontal tail buffet loads were used. In this case, 
overall agreement is improved, however, the "torsion mode peak" 
is further degraded. 

This study illustrates the importance of the coupling between 
the torsional motion and shock movement as a mechanism for pro- 
ducing severe buffeting forces. Since the wind tunnel. model was 
significantly more rigid than the airplane, its torsional motion 
could not produce a forcing function as strong as that which 
occurred in the airplane via the coupling mechanism. Also, linear 
theory unsteady aerod 3 mamics as used which do not consider the 
presence of shocks in the flow would not produce the shock motion 
induced loads « Thus, failure to predict the "torsion mode peak" 
is felt to be a result of not including the torsion-shock coupling 
mechanism in the theoretical model. 


Effect of Maneuver Transients 

The buffet prediction method was developed on the assumption 
that buffet response can be represented by a stationary random 
process at any point during a maneuver. In actuality, however, 
the maneuvers in which high intensity buffet are encountered are 
often highly transient. In such instances, predictions made at 
the limits of the maneuver should represent an upper and lower 
bound for measured flight test results. (The bounding of tran- 
sient response should not be confused with the bounds at a given 
a due to combined S3mimetric and anti S5mime trie response.) 

A pullup maneuver performed during Flight 79 for M == 0.8, 
alt. - 1490m, G.W, = 324,000N and A = 26°, (discussed in ref. 4) 
serves as an excellent example for illustrating the transient case. 
The time history of angle of attack, load factor, Mach number, and 
dynamic pressure is shown in figure 15. Flight test data analyses 
were made at a = 4.1°, 7.5°, 8.9° and 10,4°. Because of the short 
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duration of the maneuver, the data samples were only one second 
long; however, the variation of a within these samples was still 
quite large as summarized in the table below. 

VARIATION OF IN THE FLIGHT 79 
DATA ANALYSIS POINTS 


j «N0M 

a 1 

a 2 

«MAX 

” 

A a 

4.1 

2.95 

5.50 

5.50 

2.55 

7.5 

5.12 

10.10 

10.10 

4.98 

8.9 

6.20 

11.60 

11.60 

5.40 

^ 10.4 

9.00 



9.60 

12.10 

----- J 

3.10 


Predictions were made with a wing-alone representation for 
symmetric modes only at a = 6.1° and 10.18° (rigid) which are 
equivalent to a = 7.0° and 11.7° (f lexibilized) . Comparisons 
with flight test results are shown in figure 16 for the wing tip 
accelerometer. In each plot, the bounds are represented by the 
S3nnmetric predictions for a = 7.0° (lower bound) and a = 11.7° 
(upper bound), where static flexibility was accounted for. 

Flight test results are shown for a = 7.5°, 8.9° and 10.4° which 
should all fall within the bounds. Noting the range of 
variation for the three flight test points in the above table, 
it is evident that the comparisons of figure 16 verify the concept. 


CAPABILITY ASSESSMENT OF THE METHOD 


As a test of the final version of the prediction method, 
results were calculated for a wide variety of flight conditions 
and compared with flight test data. The comparisons presented in 
this section are made on the basis of RMS response and character- 
istic frequencies, No. The predictions are given in the form of 
an upper and lower bounds for all response items except pilot seat 
and C.G. accelerometers. In addition, the calculated natural mode 
shapes used in the predictions are given for all cases. 

The second volume of this report presents the predicted upper 
and lower bounds power spectra for all of the cases and response 
items given in this section. The flight test power spectra are 
shown on each prediction plot for the nominal value of angle of 
attack that most closely agrees with the flexible angle for the 
prediction. The flight test and prediction conditions are given 
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in both this voltime and Volume 2 in tabular form for all cases 
considered. The order in which the plots are given corresponds 
to that of the results given in this section. 

The conditions used for assessment of the prediction method 
are summarized as the following cases: 

Case 1: Wing alone prediction 

A = 260, M = 0.80, Alt = 6035m 

Case 2: Total airplane prediction (half horizontal tail, 

matched first wing torsion mode frequencies) 

A = 26°, M = 0.80, Alt = 6035m 

Case 3: Total airplane prediction (final method) 

A =26°, M = 0.70, Alt = 7559m 

Case 4: Total airplane prediction (final method) 

A = 50°, M = 0.85, Alt - 8383m 

Case 5: Total airplane prediction (final method) 

A = 72. 5*^, M = 0.85, Alt = 7285m 

Case 6: Wing alone prediction (final method) 

A = 50°, M = 1.20, Alt = 9053m 

Case 7: Wing alone prediction (final method) 

A = 72.5°, M = lo20, Alt = 9083m 

The methods are referred to in each case since not all cases were 
treated in the same manner. The wing alone predictions were made 
with no aerodynamic loadings on the horizontal tail or airplane 
fuselage forward and aft of the wing. In Case 1, the spanwise 
phase angles were not used whereas they were used for the wing 
alone predictions in Cases 6 and 7. In Case 2, spanwise phasing 
was not used, however, total airplane aerod 3 mamics were used. 

Case 2 was the only prediction to use shifted frequencies in order 
to match the first wing torsion modes of the full scale airplane 
and the model. Cases 3, 4 and 5 made use of the complete final 
prediction method. 

Before beginning a detailed discussion of the prediction 
results to follow, some general remarks are appropriate on the 
overall characteristics of the RMS-Nq and PSD results. The wing 
tip accelerometer results usually show the best agreement with 
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flight test whereas the horizontal tail loads are usually the 
worst. Since the horizontal tail buffet pressures were estimated 
from wing loads, this is not surprising. The predicted wing bend- 
ing moment and shear are also generally in agreement with flight 
test results. Due to the short moment arm, wing torsion is usually 
less in agreement with flight test. 

An interesting observation can be made on the c.g. and pilot 
seat accelerometers. Since they are on the centerline, upper and 
lower bounds cannot be defined for these items as they are for the 
other response items. Hence, lateral accelerations are due to 
antis 3 mimetric responses and vertical accelerations are due to 
S 3 mimetric responses. To the contrary, flight test results show 
that both vertical and lateral accelerometers each respond in both 
symmetric and antisymmetric modes. This is perhaps the strongest 
indication of the presence of as 5 mimetric modes that has been 
observed. 


Prediction Results 

The case conditions are described in tables and the prediction 
results and natural modes in figures. A cross-reference for case 
numbers, descriptive tables, natural mode figures and response 
figures is given in Table 3. 

Case 1 .- This case refers to a wing alone analysis of the 
same maneuver for which a total airplane analysis is performed 
under Case 2. The conditions were at a wing sweep of 26® and 
M = 0.80. This particular maneuver, performed during Flight 77, 
has been the subject of much of the developmental analysis per- 
formed in the previous section on the method evolution. 

The wing alone analysis made use of the chordwise pressure 
distributions without spanwise phase angle and hence does not 
completely represent the final method. The modes were calculated 
for the total airplane as shown in figures 17 for symmetric 
motions and figures 18 for antis57mmetric motions. The predicted 
RMS and Nq responses shown in figures 19 indicate reasonable 
agreement with flight test data for both the wing tip and c.g. 
accelerometers which is also backed up by the PSD results in 
Volume II (figures 1) . Wing bending moment and shear are over 
predicted in a totally unrealistic manner. Wing torsion appears 
to show excellent agreement in figure 19, which is also verified 
by the PSD results (Volume II, figure 1) with exception of response 
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near the first torsion mode frequency (about 25 Hz). 


Case 2 .- For this case, unsteady aerod 5 m.amics and estimated 
buffet loads were added for the horizontal tail in the analysis of 
the maneuver which was studied in Case 1. In addition, the scaled 
model first wing torsion frequency was matched with that of the 
airplane as indicated by the change in the RMS frequency limit 
from 38 Hz to 31 Hz. As discussed in the previous section under 
the effect of horizontal tail buffet loads, all predicted accel- 
erometer responses were improved as shown in figures 20. Since 
the modes used were the same (shown in figures 17 and 18) , the 
improvement is due only to the change in aerodynamic representa- 
tions. Wing bending moment and shear RMS are well bounded by the 
total airplane prediction. The excellent agreement with Nq is 
further emphasized by corresponding agreement for the PSD results 
(Volume II, figures 2f and 2g) . Wing torsion is under predicted 
even with the upper bounds at the intermediate angle, a = 11.7°. 
This is primarily due to the fact that the methodology cannot pre- 
dict the torsion mode peak as shown in the PSDs. The horizontal 
tail loads are over predicted as discussed previously. The PSDs 
indicate the problem is due to large inertia loads encountered 
near the horizontal tail modes at about 15, 25 and 30 Hz 
(Volume II, figures 2i, 2j and 2k). This characteristic is 
repeated in the following cases where the total airplane is 
considered. 

Case 3 .- Whereas Cases 1 and 2 were concerned with a flow 
condition in which the buffeting forces were dominated by shock 
induced loadings, this case is concerned with a lower Mach number, 
M = 0.7, where the forces are caused by a leading edge separated 
flow. The modes used are shown in figures 21 and 22 and the 
RMS-No results are shown in figures 23, The predictions are 
obtained with the final method with chordwise and s panwise phasing 
and total airplane aerodynamics. Shifting of the frequencies to 
obtain torsion mode frequency matching was not done in this or any 
of the remaining cases. 

The RMS-No results tend to follow the peaky characteristics 
of the flight test data as discussed by Benepe (ref. 2). The 
trend is further verified by the PSD results (Volume II, figure 3). 
Good agreement is shown for the accelerometer results except for 
the pilots seat. The wing loadings are well bounded with excep- 
tion of one point for wing torsion. The horizontal tail loads are 
again over predicted with large inertia loads introduced in the 
first pitch and bending modes as shown by the PSDs. (Volume II, 
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figures 3i , 3 j and 3k) . 


Case 4 .- For this case, the wing sweep was increased to 50° 
at a higher Mach number of 0,85. This case also made use of the 
final prediction method as in Case 3 above. The modes are shown 
in figures 24 and 25. 

The RMS -No predictions shown in figures 26 are higher 
relative to flight test results in that the latter tend to fall 
along the lower bounds. The characteristics seem to be well pre- 
dicted, however, as shown by the PSDs (Volume II, figures 4) in 
which the flight test data are well bounded. The horizontal tail 

loads are again plagued by the large inertia loadings due to 

horizontal tail modes (Volume II, figures 4i, 4j and 4k). 

Case 5 .- The wing sweep was further increased to 72.5° with 
about the same Mach number of 0.85. The final prediction method 
was used and the modes are shown in figures 27 and 28. 

Comparison of results for this case in figures 29 and Case 4 

shows a significant drop in response levels for both flight test 
data and predicted results. The accelerometer characteristics are 
well predicted as is more evident by the PSD results (Voltime II, 
figure 5). The wing loads are bounded for both RMS and PSD com- 
parisons. The horizontal tail loads are again over predicted 
where in this case the inertia loads in the PSDs are more pro- 
nounced near the pitch mode frequencies (about 30-34 Hz) . 

Case 6 .- For a wing sweep of 50°, the Mach number was 
increased from M == 0.85 to M = 1,2. The predictions were made 
with the final method for which wing loads only were used due to 
restrictions in the unsteady supersonic aerodynamic procedure. 

The modes used are shown in figures 30 and 31, 

The accelerometer and wing loads predictions as shown in 
figures 32 bound the flight test data very well. Comparison with 
results in Case 4 for the same sweep at a lower Mach number shows 
a significant decrease in buffet intensity. The only item which 
does not agree in either the RMS -No or PSD comparisons is the c.g. 
accelerometer. This is felt to be in part due to not having a 
total airplane simulation. 

Case 7 .- For a wing sweep of 72.5°, the Mach nimiber was also 
increased from M= 0.85 to M= 1.2. Wing alone aerod 3 mamics were 
used with the final method in the same manner as was done in 
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Case 6 above. The modes used are also given in figures 27 and 28 
(same as for Case 5). An important exception for this case was 
that the unsteady aerodynamics would not correctly handle a sub- 
sonic trailing edge configuration. Hence, the stiffness and 
damping were over predicted for this case. 

The RMS-No results in figures 33 show the effect of too much 
aerod 3 mamic damping. The responses tend to under predict as 
opposed to previously observed trends. The PSD results (Volume II, 
figures 7) further emphasize the behavior in that characteristics 
are not well predicted. The bounding of the flight test data is 
very poor for a = 8.4° and 15.5° when compared with the other 
predictions for Cases 1-6. Although this case is a failure in 
comparison with the flight test data, it illustrates the impor- 
tance of having the correct unsteady aerodynamic representation. 


Summary of the Capability Assessment 

The preceding results have illustrated the capability of the 
buffet prediction method to predict buffet response bounds on full 
scale aircraft in flight. Comparison of Cases 1 and 2 shows the 
importance of considering the total airplane when predicting the 
buffet response of an aerod 3 mamically balanced airplane in flighto 
Case 7 shows the importance of using accurate unsteady aerod 5 mamic 
forces . 

Comparison of Cases 2 and 3 shows the effect of different 
types of buffeting flow at subsonic speeds and A = 26*^. In Case 3, 
the flow at M = 0.7 is primarily of the leading edge separation 
type for which the response exhibits the peaky characteristic as 
discussed by Benepe (ref, 2). In Case 2 the flow at M = 0.80 is 
dominated by shocks on the wing which produces a far more severe 
buffeting condition. Comparison of these cases shows that Case 2 
responses are more than double those of Case 3 as indicated by 
both predictions and flight test data. 

Cases 2, 4 and 5 show the effects of increasing sweep from 
26° to 50*^ to 72 o 5° at about M - 0.80 - 0.85. Although gross 
weight and altitude are also changing, the major influence is due 
to sweep. In these cases, the decrease in buffet intensity with 
sweep as shown by flight test data, is well predicted with the 
method. 
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Case 6 shows the results for A = 50° and M = 1.2 which were 
obtained with a wing alone simulation. Comparison with Case 4 
for M = 0.85 illustrates primarily the changes which take place 
due to the Mach number changing from subsonic to supersonic. 

Again, the decrease in buffet intensity with increased Mach niimber 
is well predicted by the method. The supersonic results would be 
improved, however, if a total airplane simulation were used. 

Figure 34 shows a frequency plot of the RMS -No flight test 
data points for the wing tip accelerometer and the wing shear. 

The data are plotted according to how many points fell between 
various 10% fractional bands of the upper bounds. For example, 

14 points fell between 20% and 30% of the upper bounds for the 
wing tip accelerometer. Although the sample is small, the results 
shown in this figure establish the fact that a reasonable rela- 
tionship does exist between the predictions and flight test. The 
data points used include all of those obtained for Cases 2-6. 


CONCLUSIONS 


A method has been presented in this report for predicting 
the high intensity buffet response characteristics of highly 
maneuverable airplanes in flight. The method appears to contain 
the major ingredients responsible for airplane buffet. Since 
high intensity buffet response of aircraft is highly dependent on 
type of maneuver, atmospheric conditions, pilot characteristics 
and other items, it was felt that a peak -by-peak discussion of 
each power spectrum was not warranted. What is more important is 
that the upper bound which represents the maximum possible response 
is rarely exceeded except in special cases or for certain items. 
Likewise, the lower bound which is the minimum possible response 
level rarely falls above flight test data. These were the desired 
results . 

Predictions were presented in the form of an upper and lower 
bounds which were verified by extensive comparison with flight 
test data. The changing spectral characteristics of flight test 
data were shown to be predictable with the method. Static aero- 
elastic effects were shown to have a significant impact on the 
predictions. The horizontal tail loads ^ere found to have a 
significant influence on wing and fuselage responses, hence, the 
total airplane should be considered. The means for estimating the 
tail buffet loads from wing data seems to be correct in concept. 
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In constructing the buffet forcing function, it was shown 
that the chordwise and s panwise phasing in the pressure data was 
important. The influence of wing torsional motion in the existing 
prediction model was found to be insignificant. This was felt to 
be due to not being able to consider normal shock motions in the 
development of the unsteady aerodynamic forces due to structural 
response. It was also shown that maneuver transient effects could 
be quite pronounced but that the response could be bounded by 
considering the extremities of the angle of attack range achieved 
during the maneuver. 

Although the prediction of buffet response as described in 
this report is a formidable task, it is felt that the mechanics 
of high intensity buffet are much better understood as a result 
of the study. Since the buffet forcing function was the most 
difficult input data to obtain, the next step would be to develop 
a means for calculating this data through a theoretical approach. 
The fact that linear theory can be used to estimate the horizontal 
tail loads due to buffeting is strong evidence that such a develop' 
ment is possible. 
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TABLE 1.- BY7 ANALYSIS/TEST COMPARISON OF 

F-lllA SYMMETRIC VIBRATION FREQUENCIES, 
A = 26 . 0 ® 


GW=201,000 N =45,180 lb GW=211,500 N =47,500 lb 



Analysis 

Results 

! Test 

Results 

(A/P #12) 

Order 

# Freq . 

Mode 

' Order // 

Freq. 

Mode 

! 

1 

1 

5.17 

1st Wing Bend 

^ 1 

5.2 

! 

1st Wing Bend \ 

1 




' 2 

7.9 

j 

Longitudinal 
Wing Bending ; 

2 

8.74 

1st Fuse Bend 

1 3 

8.6 

1st Fuse Bend 

1 

1 3 

15.08 

H. Tail Bend 


13.6 

H. Tail Bend i 

i 

j 

4 

16.12 

2nd Wing Bend 

5 

16.9 

2nd Wing Bend ! 

■; 5 

20.98 

2nd Fuse Bend 



J 

1 

6 

25.48 

Wing Torsion 

6 

25.2 

Wing Torsion ' 

i 

; 7 

28.45 


1 



! 8 

31.48 


1 

1 



9 

33.37 

H. Tail Torsion 

7 

34.4 

H. Tail Torsion 

10 

35.64 
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TABLE 2.- BY7 ANALYSIS/TEST COMPARISON OF F-lllA 
ANTISYMMETRIC VIBRA.TION FREQUENCIES, 

A = 26. 0® 


GW=201,000N =45,180 lb GW=211,500N =47,500 lb 



Analysis 

Results 

Test 

Results (A/P #12) 

Order 

# Freq. 

7.56 

Mode 

Order # 

Freq . 

Mode [ 

- * ■ 

1 

1st Wing Bend 

1 

7.6 

i 

1st Wing Bend | 




2 

9.3 

1 

Longitudinal 
Wing Bending 

2 

9.93 

Lat Fuse Bend 

3 

9.9 

Lat Fuse Bend 

3 

11.54 

V. Fin Bend 

1 

11.7 

This mode omitted 
(Freq given is 
for different GW) ' 

4 

12.52 

H. Tail Bend 

4 

13.3 

H. Tail Bend j 



(with V. Fin 
Interaction) 

i 

1 5 

t 

16.2 


j 

: 5 

18.61 

2nd Wing Bend 

f 

: 6 

17.5 

2nd Wing Bend 

6 

20.79 

2nd Lat Fuse Bend 




7 

23.25 





8 

24.68 

Wing Torsion 

■ ^ 

25.4 

Wing Torsion 

9 

25.41 

Rudder 

! 7 

1 

23.6 

Rudder 

10 

25.92 

V. Fin Torsion 

i 9 

28.0 

1 

V. Fin Torsion j 

11 

28.26 

Tail Interaction 

; 10 

29.2 

V, Fin Torsion i 



plus 2nd Wing Bend 

i 


plus 2nd Wing Bend 

12 

31.68 

H. Tail Torsion 

i 

i 11 

31.0 

H. Tail Torsion i 



plus V. Fin 

1 


plus V. Fin j 

1 


Torsion 

j 

I 12 

32.7 

Torsion ' 

1 



! 
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TABLE 2 . - Continued 


GW=201,000 N =45,180 lb GW=211,500 N =47,500 lb 



Analysis 

Results 

Test 

Results (A/P #12) 

i Order 

# Freq . 

Mode 

Order # 

Freq. 

Mode 

1 

; 

i 

32.79 

3rd Lat Fuse Bend 




^ 14 

i 

37.09 

H. Tail Torsion 

13 

37.3 

H. Tail Torsion 

i 

i 15 

1 

i 

1 

39.02 

Higher H. Tail Mode 

14 

43.8 

Rudder Bending 
and Torsion 

1 



15 

45.0 

Rudder Bending 
and Torsion 
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TABLE 3.- CROSS REFERENCE FOR CASE 

NUMBERS AND TABLES AND FIGURES 


CASE 

DESCRIPTIVE 

TABLE 

NATURAL 
MODE FIGURES 

! RESPONSE i 
FIGURE ‘ 

1 

4 

17(sym), 18 (anti) 

19 

2 

5 

17(syTn), 18 (anti) 

20 

3 

6 

21(syin), 22 (anti) 

23 

4 

: 7 

24(syni), 25 (anti) 

26 

5 

; 8 

27 (sym) , 28 (anti) 

29 

6 

9 

30(sym), 31 (anti) 

32 

7 

10 

27 (sym), 28 (anti) 

33 


TABLE 4.- CASE 1, WING ALONE, A = 26°, M = 0.80, 
ALT = 6035m (19,800 ft) 


CASE 1 FLIGHT 77 RUN SC-R 



FLIGHT 

TEST CONDITIONS 


PREDICTION 

CONDITIONS 

Wing Sweep 

25.6° 




26° 


Mach 


0.80 - 

D.78 



0.80 


Altitude 

6035m (19,800 ft) 



6035m (19,800 ft) 

Gross Weight 

266,004N (59,800 lb) 


266,004N (59,800 lb) 

1 POINTS ANALYSED 



AT 

“1 

“2 

“MAX 

“NOM 


* 

“FLEX 

“RIG 

2 

4.22° 

5.98° 


5.1° 




2 

6.80° 

7.12 




6.6° 

6.1° 

2 

2 

8.15„ 

10.35^ 

9.35„ 

12.90. 

• 

8-9o 



10.18° 

2 

12.70° 

14.65° 

14.95° 



14.4 

13.26° 

3 

11.05° 

14.95 

14.95 

13.0 


" 



* ^^^FLEX _ ^RIG 
^^«*RIG “FLEX 


0*920 as obtained from Figure 4 



TABLE 5.- CASE 2, TOTAL AIRPLANE (HALF HORIZONTAL TAIL) 
A = 26°, M = 0.80, ALT = 6035m (19,800 ft) 


CASE 2 


FLIGHT 77 RUN SC-R 



FLIGHT TEST CONDITIONS 

PREDICTION CONDITIONS 

Wing Sweep 
Hach 
Altitude 
Gross Weight 

25.6° 

0.80 - 0.78 
6035m (19.800 ft) 
266.004N (59,800 lb) 

26° 

0.80 

6035m (19,800 ft) 
266.044N (59,800 lb) 

1 POINTS ANALYSED 


A T 

“1 

«2 

O' MAX 

“NOM 

* 

“FLEX 

“RIG 

2 

2 

2 

2 

2 

3 

4.22^ 

6.80^ 

8.15° 

10.35° 

12.70° 

11.05 

5.98° 

7.12^ 

12.90° 

14.65° 

14.95 

14.95° 

14.95 

13.0 

6.6° 

11.1° 

14.4 

6.1° 

10.18° 

13.26 


* ^^“FLEX 

Cl a RIG 


°RIG 

“flex 


= 0.920 as obtained from Figure 4 


TABLE 6.- 


CASE 3 , TOTAL AIRPLANE (FINAL METHOD) , 

A - 26°, M = 0.70, ALT = 7559m (24,800 ft) 


CASE 3 FLIGHT 48 RUN 6 



FLIGHT TEST CONDITIONS 

PREDICTION CONDITIONS 

Wing Sweep 
Mach 
Altitude 
Gross Weight 

26.6° 

0.70 - 0.68 
7559m (24,800 ft) 
294.472N (66,200 lb) 

26° 

0.70 

7559m (24,800 ft) 
293,138N (65,900,1b) 

POINTS ANALYSED 


AT 

“1 

«2 

"MAX 

a 

* 

“FLEX 

“RIG 

1 

1 

1 

1 

1 

8.72° 

9.70^ 

10.30^ 

11.15° 

14.25 

9.55° 

10.75^ 

11.75° 

13.55^ 

16.60 

- 

8.8° 

9-So 

10. 7„ 
11.8° 
14.6 

10.7^ 
11.8„ 
12. 8„ 
17.1 

10.2^ 

12.2^ 

16.3 


* ^ “RIG 

^LOrig 


0.950 as obtained from Figure 4 


FLEX 




TABLE 7 . - CASE 4 , TOTAL AIRPLANE (FINAL METHOD) , 

A = 50°. M = 0.85, ALT = 8382m (27,500 ft) 

CASE 4 FLIGHT 61 RUN R227 



_ — -- — — n 

FLIGHT TEST CONDITIONS j 

PREDICTION 

CONDITIONS 

Wing Sweep 

49.1° 



50° 


Mach 


0.82 - 

0.79 


0.85 


Altitude 

8382m (27,500 ft) 


8382m (27,500 ft) 

Gross Weight 

330,948N (74,400 lb) 

331,392N (74,515 lb) 

1 POINTS ANALYSED 







* 


AT 

“1 

«2 

“MAX 

“NON 

“FLEX 

“RIG 

1 

7.10° 

9.25° 


7.9° 


_ 

1 

8.05° 

10.10° 

- 

8.9° 

8.9° 

8.1° 

1 

10.10° 

10.80° 

- 

10.0° 

- 

- 

1 

10.60° 

12.70° 

- 

11.1° 

11.1° 

10.2°' 

1 

12.90° 

14.60° 

14.60° 

13.1° 

14.4° 

13.2° 


* ^ ^ FLEX „ _^RIG _ 0.920 as obtained from Figure 4 

^LotRIG “flex 


TABLE 8 . - CASE 5 , TOTAL AIRPLANE (FINAL METHOD) , 

A = 72.5°, M = 0.85, ALT = 7285m (23,900 ft) 


CASE 5 FLIGHT 48 RUN 7-R 



FLIGHT 

TEST CONDITIONS 


PREDICTION 

CONDITIONS 

Wing Sweep 

72.2° 




72.5° 


Mach 


0.89 - 

0.84 



0.85 


Altitude 

7559m (24,800 ft) 



7285m (23,900 ft) 

Gross Weight 

265.559N (59.700 

lb) 


268.673N (60,500 lb) 

1 POINTS ANALYSED 









ic 


AT 

“1 

“2 

“MAX 

“NOM 


“FLEX 

“RIG 

1 

7.15° 

8.65° 


7.8° 


7.8° 

7.1° 

1 

8.65° 

10.00° 

- 

9.4° 


- 

- 

1 

10.75° 

12.20° 

- 

11.1° 


11.1° 

10.2° 

1 

14.15° 

16.15° 

- 

14.4° 


14.4° 

13.3° 

1 

17.90° 

18.90° 

19.35° 

17.7° 





* ^^“FLEX „ “rig 
CLa^jg “flex 


0.890 as obtained from Figure 4 




TABLE 9 . - CASE 6 , WING ALONE (FINAL METHOD) , 

A = 50°, M = 1.20, ALT = 9053m (29,700 ft) 


CASE 6 FLIGHT 48 RUN 4 



FLIGHT 

TEST CONDITIONS 


PREDICTION 

CONDITIONS 

Wing Sweep 

49.8° 




50° 


1 Kach 


1.20 - 

1.15 



1.20 


Altitude 

9053m (29,700 ft) 



9053m (29,700 ft) 

Gross Weight 

261,111N (58,700 

lb) 


261,778N (58,900 lb) 

1 POINTS ANALYSED 









* 


AT 

"1 

ao 

"MAX 

"NOM 


"FLEX 

"rig 

1 

4.70° 

5.50° 


4.9° 




1 

8.20° 

9.80° 

- 

8.6° 


- 


1 

12.10° 

13.70° 

~ 

12.4° 


12.4° 

10.2° 

1 

13.70° 

13.90° 

15.0° 

13.7 


16.1° 

13.2° 


* = 0.823 as obtained from Figure 4 

"flex 


TABLE 10.- CASE 7, WING ALONE (FINAL METHOD), 

A = 72.5°, M = 1.20, ALT = 9083m (29,800 ft) 


CASE 7 FLIGHT 48 RUN 5 



FLIGHT 

TEST CONDITIONS 

■ 

PREDICTION 

CONDITIONS 

Wing Sweep 

72.2° 




72.5° 


Kach 


1.20 - 

1.16 



1.20 


Altitude 

9083m (29,800 ft) 



9083m (29,800 ft) 

Gross Weight 

274,455N (61,700 lb) 


268,673N (60,500 lb) ’ 

POINTS ANALYSED 









* 


AT 

"1 

"2 

"MAX 

"NOM 


"FLEX 

"rig 

1 

4.80° 

4.80° 


4.8° 




1 

8.00° 

8.80° 

- 

8.1° 


8.1° 

7.1° 

1 

11.30° 

12.70° 

- 

11.6° 


11.6° 

10.2° 

1 

14.95° 

16.75° 


15.1° 


15.1° 

13.4° 


* ^^"FLEX „ <^RIG 
Cl a RIG "FLEX 


0,837 as obtained from Figure 4 






Figure 1.- Buffet prediction method flow diagram 
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(a) Relationship between the NASA ARC 
spectral density analysis and the 
pressure transducer locations 
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(b) Spline fit of model pressures for 

M=0.80, a=10.18° and f , ,=48.0 Hz 

model 


Figure 5.- Conversion of wind tunnel PSD's 

into complex pressure distributions 




a Cl 

RIGID PREDICTION FLEXIBLE PREDICTION 


Figure 6.- Rigid and flexible wing bending for 
A = 26°, M = 0.8, alt = 6035m 




POWER SPECTRAL DENSITY C6D»*2/HZ 
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FREQUENCY C HZ ) 


Figure 7.- Predicted bounds on wing tip accelerometer PSD 
for A = 26°, M = 0.8, alt = 6035m 
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POVER SPECTR^L DENSITY CM-N5*»2/HZ 


4 > 




FREQUENCY t HZ ) 


Figure 8.- Predicted bounds on wing bending PSD for 
A = 26°, M = 0.8, alt = 6035m 




O PREDICTED UPPER BOUNDS 
PREDICTED LOWER BOUNDS 
FLIGHT TEST RESULTS 


Figure 9 .- Effect of horizontal tail loads on wing 
shear for A = 26°, M - 0.8, alt = 6035m 





o predicted upper bounds 

B PREDICTED LOWER BOUNDS 

A O FLIGHT TEST RESULTS 


Figure 10.- Effect of horizontal tail loads on wing tip 

accelerometer for A = 26®, M = 0.8, alt = 6035m 
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o PREDICTED UPPER BOUNDS 
□ PREDICTED LOWER BOUND 
A FLIGHT TEST RESULTS 


Figure 11. -- Effect of horizontal tail loads on c.g. 

accelerometer for A = 26®, K - 0.8, alt = 6035m 


57 



POWER SPECTRAL DENSITY (C**2)/H2 POWER SPECTRAL DENSITY (C**2)/HZ 
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FREQUENCY (RE) 



5 10 15 20 25 30 33 
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NO PHASING 


SPANtfISE PHASING 



CHORDWISE PHASING 



Figure 12.- Effect of buffet |?ressure spanwise and chordwise 
phasing on wing alone response PSDs for 
A = 260, = 0.8, alt = 6035m 

(a) ct = 6,6^ 
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© CHORDWISE PHASE ANGLE ONLY 
A FLIGHT TEST DATA 
E? NO PHASE ANGLE (SPANWISE LOADING ONLY) 
6 SPANUISE PHASE ANGLE 
^ CHORDWISE AND SPANWISE PHASE ANGLES 


Figure 13.- Effect of buffet pressure spanwise and chordwise 
phasing on wing alone response RMS and Nq results 
for A = 26°, M == 0.80, alt = 6035m 
(a) accelerometer response 
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Figure 13.- (b) wing loading response 
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Figure 19*- Response predictions for case 1 
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Figure 19.- (continued) 
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Figure 1 9 -- (continued) 
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Figure 20. -(continued) 
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Pigure 20,- ( continued j 
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Figure 23.- Response predictions for case 3 
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figure 2 3. -(continued) 
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Fipi;ure 24. -(continued) 






98 


A=50.0° and G.W.=331 ,392N 


Mode 4 12*290 Hz 


Mode 5 


Mode 6 



t 




Figure 25 (continued) 
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Figure 26.- Response predictions for case 4 
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Figure 26 .-(continued) 
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Figure 28 .-(continued) 
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Figure 29.- Response predictions for case 5 
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Figure 29 .-(continued) 
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Figure 29. -(continued) 
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Figure 29 . -(continued) 
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Figure 52.- Response predictions for case 6 
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Pigure 32. -(continued) 
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Figure 32*- (continued) 
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Figure 55.- Response predictions for case 7 
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Figure 35.- (continued) 
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FLT 48 RUN 5 
M - 1.2 
h - 9083 M 
GW - 268,673 M 
SWEEP - 72,5° 
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Figure 33 .-(continued) 





